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We investigate far-infrared properties of strain engineered vanadium dioxide nanosheets through 
epitaxial growth on a (100) j?, Ti02 substrate. The nanosheets exhibit large uniaxial strain leading to 
highly uniform and oriented cracks along the rutile c-axis. Dramatic anisotropy arises for both the 
metal-insulator transition temperature, which is different from the structural transition temperature 
along the cr axis, and the metallic state conductivity. Detailed analysis reveals a Mott-Hubbard 
like behavior along the rutile cr axis. 
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I. INTRODUCTION 

During the past five decades, vanadates have been intensely investigated from the point of view of clarifying 
the physics of metal- insulator transitions (MIT). Vanadium dioxide (VO2) is particularly intriguing, exhibiting a 
MIT with a conductivity decrease of over five orders of magnitude^. This is accompanied by a transition from a 
high temperature metallic rutile structure to an insulating monoclinic structure below the transition temperature. 
In bulk, the structural transition temperature, 15""^, equals the metal-insulator transition temperature, Tj/j'^, i.e. 
rj^buik ^ rj^buik ^ 340^^2 However, VO2 is not yet fully understood regarding the Mott-Hubbard or Peierls like nature 
of its MIT, though i t is now generally accepted that a complete description requires explicit consideration of electron- 
electron correlation^^. This is motivated in particular by the observed difference between the MIT and the structural 
transition temperatures, which provides significant evidence that the correlated conductivity behavior is independent 
of the structural phase^^. Additionally, recent experimental studies highlight the multifunctional potential of VO2 
including current or photoinduced phase changes (with the commensurate large tuning of the dielectric function) as well 
as integration with other technologies for applications including light detection and memory-based metamaterial^^'l^. 

Strain engineering enables an additional degree of control of technologically relevant properties and provides a 
discriminatory capability towards obtaining fundamental insight into the microscopic origin of the macroscopic char- 
acteristics. Strain ha s been used to modify the MIT temperature in VO2 through direct applicatiorl^ or substrate- 
dependent applicatiorP^EU of stress. For the present measurements, we utilize highly strained epitaxial (lOO)j^ VO2 
thin films. The far-infrared conductivity is measured using non-contact polarization sensitive terahertz time-domain 
spectroscopy (THz-TDS) providing direct access to the coherent quasiparticle response along b^ or cr by simply 
changing the sample orientation. Subsequent analysis indicates that the strain engineered tuning of vanadium aig 
and Bg orbitals controls the MIT transition temperature along c^, which is different from the structural transition 
temperature. The origin of the observed conductivity anisotropy is also discussed, along with the potential for a 
technological application of strain engineered VO2 thin films as temperature switched far-infrared polarizers. 

II. EXPERIMENTS 

The ~ lOOnm and ~ 250nm thick VO2 films we investigated were deposited on rutile (lOO)ii; Ti02 substrates by 
temperature optimized sputtering from a vanadium target, using the reactive bias target ion beam deposition technique 
in an Ar-|-02 gas mixture; details of the growth conditions can be found elsewhere^. The samples morphology was 
characterized by optical, atomic force (AFM) and scanning electron microscopies (SEM), and the film microsctructure 
was analyzed by temperature dependent X-ray diffraction (XRD). 

THz-TDS is a non-contact method to measure far-infrared conductivity. The conductivity anisotropy is easily 
determined from transmission of THz pulses for different sample orientations. In our case the output of a IkHz 35fs 
Ti:Sapphire amplifier is used to generate nearly single-cycle THz pulses via optical rectification in a ZnTe crystal. We 
then employ a standard THz Time Domain Spectroscopy (THz TDS) setup to measure the transmission in the THz 
frequency range ~ 0.2 — 2.0THJ1^. By changing the sample orientation with respect to the THz pulse polarization 
we are able to monitor the transmission along different crystal axes. 

(lOO)i^ Ti02 has a large refractive index anisotropy in the THz range, partly caused by its characteristic rutile 
structure. Such an anisotropy in the (100)^ Ti02 substrate makes it easy to distinguish transmitted THz signals 
with polarization parallel to c^ from THz signals with polarization perpendicular to b/j (Fig. [l|a)). Temperature 
dependent experimental characterization of the Ti02 substrate in the THz range, using our THz TDS setup, led to 
values of 

He ~ 12.4 - 6 X lO^^AT + 1.2i 
Tio ~ 9.1 - 3 X 10-*AT + 0Ai 

for the refractive index along c_r and b^, respectively, in accordance with previous result J^^. Here, AT stands for 
the temperature deviation with respect to room temperature. This approximation to the Ti02 complex refractive 
index holds between room temperature and 400K; it does not take into account the temperature dependence of the 
imaginary part, which is negligible for our purposes. 

The conductivity of the VO2 film is extracted using the Fresnel equations, after experimental determination of the 
ratio of the THz transmission of the film to the THz transmission of a bare Ti02 substrate, used as a reference. 
THz TDS is thus a non-contact conductivity measurement, which allows one to quickly characterize the conductivity 
anisotropy in many samples with different thicknesses and substrates. 

III. RESULTS 

A. Characterization of the VO2 nanosheets 
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Figure [2]ja) is an AFM phase image of the 250nm film, showing periodic buckhng and cracking of the film parallel 
to c/j, with inter-crack spacings on the order of 1/im. The depth profile reveals nanosized ridges, ~ 15nm high, near 
the cracks. The SEM image in Fig. [2jb) confirms the ~ 250nm thickness of the film. It shows that the cracks' depth 
matches the film thickness and it gives an estimate of about SOnm for their width. Such nanocracks were not detected 
by the AFM due to lack of tip sensitivity. 

The optical images of the lOOnm (Fig. [sf^a)) and the 250nm (Fig. [sf^b)) films confirm the ^ Ifim period spacing 
of the cracks along Cj^. Also, the observed uniformity in the distribution of cracks is an indication that our films are 
strained in a highly homogeneous and oriented fashion. This crack distribution enables a comparison of our results 
with those of VO2 nanobeams, aligned along c jSEIHlol^ although the dimensionality of our nanosheets gives access 
to the properties of strained VO2 along more than one axis. Such cracking uniformity was achieved by optimizing 
the growth conditions of the samples, in particular the growth temperature. In the case of our (100)^ VO2 films the 
growth temperature was set to ~ 500°C. Films grown at ^ 450°C show cracking along different directions, which 
makes the characterization of the films along cr less straightforward. All subsequent analysis is analogous for both 
the lOOnm and the 250nm thick films. Unless otherwise specified the results will refer to the lOOnm thick sample. 

Detailed room temperature XRD of the VO2 film confirms the single crystal nature of the sample and yields the 
lattice parameters of the material (Fig. |4][a)). Comparing these values with those in the rutile phase of bulk V02^ 
yields mismatches of —0.83% along a^j, —2.17% along hu, and 1.41% along cu, indicative of a compressive strain along 
a,ji and b^, as opposed to a tensile strain along cr. Such strain values along c^ are comparable to those achieved 
in one-dimensional VO2 nanobeamJi^. Bulk rutile Ti02 has larger lattice constants than bulk rutile VO2 along all 
directions so one would expect that both ha and c^, would expand in (lOO)ii VO2 films grown on a {100)r Ti02 
substrate. However, our nanosheets show that the expansion along cr surpasses the substrate clamping effect due to 
the bfl-axis, leading to a compressive strain in the ab ^-plane. This behavior is also observed in the 250nrn thick film, 
with mismatches of -0.68% along a^sO"™, -1.94% along b|f""™ and 0.86% along c|5°"". 

Temperature dependent XRD results (Fig. [3](c)) show that a small a^j-axis expansion, ~ 0.1%, occurs during the 
film's structural transition. This transition occurs at a temperature T/f" - 340^ - T^f'', the same as in bulk, 
and shows the expected hysteric behavior. (From here on Tgt will refer to both bulk and film structural transition 
temperatures.) The ag-axis expansion in our strained sheet can be compared to that in bulk VO2, estimated as 



"^fcjv/ ~ ^'^ ~ 0.8%P22l^ where bM is the equivalent of in the bulk monoclinic insulating structure. The order of 
magnitude difference between the lattice parameter variation in the film and that in bulk stems from the clamping 
effect of the rutile Ti02 substrate. This is consistent with the large strain observed in room temperature XRD, caused 
by the strong substrate clamping effect. The structural change exhibited at 3iOK by our strained VO2 nanosheets is 
thus smaller than the change observed in bulk V02. Since VO2 shares the Ti02 rutile structure at high temperatures, 
the reduced structural change is expected to have an impact mainly on the low temperature VO2 film's structure. 
The temperature dependent (200)_r 29 plots from which the data in Fig. 3]jc) were extracted are presented in Fig. 
^h). These raw data provide further support to our observation that the c' 



amping effect due to the substrate is very 



strong, thereby preventing the development of significantly different structural phases in the strained VO2 film. 

B. THz Time Domain Spectroscopy 

As shown in Fig. [ij upon increasing the temperature from the insulating to the metallic phase the THz peak 
transmission in the lOOnm (100)^ VO2 film drops by ~70% along the c^^-axis and by ^15% along the b/j-axis, 
and the THz peak transmission in the 250nm (100)^ VO2 film drops by ^85% along the C/j-axis and by ^--^15% 
along the b^j-axis. Figure [ijb) shows the transmission anisotropy in our 250nm VO2 sample, normalized to its low 
temperature value along each axis (bij and cr). This representation highlights the dramatic difference between the 
low temperature transmission along b^ and that along c^, thereby illustrating the potential of strained VO2 films as 
temperature tunable THz polarizing beamsplitters. 

Figure [sjja) shows the temperature dependent far infrared conductivity, obtained from the transmission data, for 
the lOOnm (100);? VO2 nanosheet. The conductivity along cr shows a clear transition from the insulating to the 
metallic state with a narrow hysteresis, which is indicative of the high quality of the film. In the metallic state 
the conductivity is acj^ ''^ 5650(J7cto)~^, comparable to bulk single crystal valued. The MIT along c^; occurs at 
a temperature T^f/ — 365° -fC. T^^j is significantly larger than both the structural transition temperature and the 
bulk MIT temperature, Tst ~ T^n'' ~ 3A0K. Our VO2 films therefore exhibit, along Cr, a, ^ 25K temperature 
difference between the structural and the metal-insulator transition temperatures. The combination of the quasi 
three dimensionality of our nanosheets, which enables a direct measurement of the strain along the three crystal axes 
through XRD analysis, with the polarization sensitivity of THz spectroscopy is the key to identifying this distinction 
between the two transition temperatures. 

The conductivity along b^ also exhibits a transition (see inset of Fig. S^a)), which occurs at T'^f^j ~ 3A0°K. Along 
this direction we therefore observe that the structural and metal-insulator transition temperatures are the same, 
^^Mi ^ Tst ~ Tt/f . However, the conductivity along hR reaches a high temperature value about 30 times smaller 
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than the high temperature conductivity along cr- This strong conductivity anisotropy will be addressed later in the 
text. 

THz TDS results for the 250nm thick sample (Fig. |6]) indicate that its transport properties are similar to that of 
the lOOnm sample. In particular, the high temperature conductivity along c/j remains as good as in bulk VO2 single 
crystalf'^ while ~ 365K > T'^^f^ ■ The high temperature conductivity along is very low, < 100(57cto)~-'^, and 
the transition temperature can only be estimated at ~ iAQK, consistent with T^f^ and Tst and in line with 
what is observed in the lOOnm sample (Fig. [sf^a)). 

IV. DISCUSSION 

Understanding the contribution of the V3d orbitals to the electronic properties is crucial in order to explain the large 
material anisotropy in Tm A splitting of the 5- fold degenerate 3d states occurs due to the octahedral coordination 
of the V atoms, resulting in a higher energy doubly degenerate Cg level and a lower energy triply degenerate t2g 
level. Trigonal distortion further splits the t2g levels leading to an upshifted doubly degenerate state, responsible 
for conduction in the ab^j-plane, while downshifting a non- degenerate cj^-oriented aig state (Fig. [5|b)). Recent 
cluster Dynamical Mean Field Theory (cDMFT) calculation^, which include the effect of a ±2% strain along Cij, 
have demonstrated that a tensile strain along narrows the aig derived band and leads to a compressive strain in 
the ab/j-plane, which uplifts the band (Fig. [sjc)). In the Mott picture the energy increase of e" electrons reduces 
the screening of electrons residing in the aig band, thereby enhancing the effect of correlations (i.e. increasing the 
screened Hubbard U). This results in an increase of the insulating band gap which opens, below Tj\,f/, between the 
bonding aig and the anti-bonding levels, therefore leading to an increased Tmi along cpi^^^. 

The epitaxial strain in our film can be decomposed into a uniaxial tensile strain along cr and a uniaxial compressive 
strain along hR. Along c^, T^ff - 365fs: > T^J'f , in hne with previous experimental resultJ^, while CTch remains 
comparable to the best single crystal valueJ^I. These results agree with the Mott picture above, where the increased 
lattice spacing along c^ further increases the ratio of Coulomb repulsion to the inter-atomic hopping integral, thereby 
destabilizing the metallic phase and increasing T^fj. The Peierls picture predicts a decrease in T'^J^j with tensile strain 
along Cjj, thus failing to describe our results along that axi^. Also contrary to the Peierls-driven MIT scenario along 
cr is the fact that the structural transition occurs at a temperature Tst ^ 340if, 2^K lower than T'^fj'- in a Peierls 
picture both transitions would be expected to occur at the same temperature. In contrast to what happens along c^, 
the results along b/j, where = Tst, are compatible with a Peierls-driven MIT scenario. 

As mentioned above, o ur sam ples are cracked along c^. The occurrence of such cracks is common in VO2, in both 
bulk and strained sample^^- and this has prevented accurate measurements of the DC conductivity in this material 
in the direction perpendicular to Cr. Hindered quasiparticle motion along hR due to the presence of Cr oriented 
cracks is the most likely explanation for the reduced value of ah^. in our films. However, a different scenario could 
arise in the THz range. Due to the small value of the far-infrared carrier mean free path (~ A-®) compared to DC, 
the ^ Ips long THz field should be able to couple to the material along hR despite the cracks along cr. According to 
this picture, the uniaxial compressive strain along b^ would play a significant role in reducing crj^, compared to cTc^, 
while keeping T^fj = T^fj'' — Tst- Given the orbitals orientation, conductivity in the abfl-plane is mediated by the 
Cg orbitals whose energy is controlled by the overlap between 02p and V3d orbitals'^''. Low conductivity behavior 
in the compressed ab/j-plane for T > Tmi would then arise from the states being higher in energy than in the 
unstrained case, which reduces their overlap with the Fermi level (Fig. pp). VO2 nanosheets that remain uncracked 
while maintaining a high level of strain along hR must be investigated in order to clarify the role of the orbitals 
on the high temperature value of (Tf,^, thereby also clarifying the nature of the MIT along that axis. 

Our strained VO2 films on a Ti02 substrate have been shown to enable the separation of the far-infrared signal's 
polarization components both in time and intensity. The temperature dependent THz conductivity of strained VO2 
nanosheets (Figs, [sjja) and |6| shows that incident light polarized parallel to b^ is transmitted through VO2 at 
both low and high temperatures, whereas incident light polarized parallel to cr is transmitted through VO2 at low 
temperatures but reflected at high temperatures. Our strained VO2 single crystal nanosheets can thus be thought of 
as temperature switchable far-infrared polarizing beamsplitters. 

V. PHASE DIAGRAM 

Finally, it is interesting to map the strain of our film on to the VO2 phase diagram in Fig. [7j constructed based on 
past studies of VO2 as a function of temperature, pressur^i^MI and dopin^2nEI, in general, it is not trivial to map 
the effect of doping onto that of pressure. This was initiated by Pou get et al., who found a one-to-one correspondence 
between Cr-doping and the application of uniaxial stress along [llOjflBMSl^ Other doping experiments tend to indicate 
that doping ions with radii smaller than V^^ lead to an increase in Tmi whereas those with larger radii have the 
opposite effect This is similar to what happens in V2O3, a canonical Mott insulator '^^^ 
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Figure [t] compiles previous experimental results obtained on VO2 under different pressur j' ^^ l ^^ l anddopin^S^Ml 
conditions. P = corresponds to atmospheric pressure. The application of hydrostatic pressur^^^^oEI] cannot be 
easily interpreted in terms of its influence on the conductivity nor on Tmi along individual axes of the crystalline 
structure, which corresponds to the main focus of our study. We therefore do not include hydrostatic pressure data in 
the phase diagram of Fig. [7] The dashed white lines delimit the region of the phase diagram where uniaxial pressure 
was applied along [001]_r (up to 12kbar), both for compression (P < 0) and tension (P > 0), and along [HO];?, (up 
to ~ 1.2kbar), only for compression (P > 0). Within the rectangle, the dashed red line represents T'^fj versus stress 
along cu, [OOlJii;, while the solid red line represents T^fj versus stress in the abij-plane, [llOj^j. The dashed black 
lines separate different insulating phases within the low temperature monoclinic phase. The dash-doted line between 
phases M2 and M4 is a conjecturd^oisai^ Note that apart from the dashed red line, relative to stress along [001]^, all 
the information in the diagram refers to stress along [110]/^. The white oval indicates the phase space characteristic of 
our (100)h VO2 samples. The uniaxial stress on our films along [110]^ has been estimated, respectively for the lOOnm 
and the 250nm samples, to ~ 20.9kbar and ~ IS.bkbar, based on the measured —2.32% and —2.06% mismatch along 
(llO)fl, and on a Young modulus of ^ 900fc6a7''*^. These uniaxial stress values are higher than ever reported, to our 
knowledge, for uniaxial stress on VO2. 

The dependence of Tmi on doping, studied in Yi-x^x02 compounds, can be approximately mapped onto its 
dependence on applied pressure. Reducti on of V^"^ is achieved using M = Nb, Mo, W, Ta, Re, Ir, F, Ti, Os, Ru, Tc, 
etc., with formal charges of +4, +5 or +^ 23 | 32 | 34 | 43| t54] rpj^g ^g^^^ 

on Tmi is similar to that of negative (compressive) 
stress along [001]r, yielding dTMi/dx = -0.3 ~ - 28K/at. %M Oxidation of is achieved using M = Cr, Al, Fe, 
Ga, Ge, Sn, Mn, Co, etc., with formal charges of +3 or +42213ilt34139143i55HlII. The effect on Tmi is similar to that of 
positive (tensile) stress along [001]^, yielding dTMi/dx = 13K/at.%M. 

Our results, indicated by the white oval in the phase diagram of Fig. [Tj are comparable to those of Everhart et ai, 
where an anisotropy of about two orders of magnitude is observed in the metallic phase of bulk VO2 single crystals 
doped with iron at 0.076%P^. Our estimate of the stress puts our samples in a region of the phase diagram where Tmi 
increases with stress, and where a two phase behavior is expected above Tmi- This is consistent with our observation 
of both insulating (T < P^fj) and metallic (P > T^fj) values of ctc^ in the same rutile structure, above Tst- The 
independence of the conductivity behavior from the structural phase has been reported befor^^lE]^ though relative to 
monoclinic structures. As for the low temperature structure, M4 or rutile seem to be the most likely can didates, 
although we cannot distinguish between the two. Ml, M2 and M3/T are all structurally too far from rutil^SSlMlMJ 
to lead to the results shown in Fig. [3p, and they are located in different regions of the phase diagram, as seen in 
Fig. [Tj These observations, along with the strong conductivity anisotropy exhibited by our strained nanosheets, are 
indicative of a more complex behavior of VO2, beyond the currently accepted doping and strain dependence of its 
structural and transport properties (Fig. [?]). 

Goming back to the conductivity measured in our strained nanosheets, it is important to point out that its 
anisotropy i s unexp ectedly large compared to previous experimental observations and theoretical calculations on 
VO2 s ample^i^'^Si'Sll xhe anisotropy in the DC conductivity, or , in undoped VO2 samples is generally > 1 for 

T < Tmi but can take many different values for T > Tmi (Table M. Differences in sample quality and stoichiometry 
as well as in conductivity measurement techniques surely affect the results, but variations in geometry and internal 
strain / cracking are likely to also have an effect on the anisotropy of metallic VO2. 

Tables [I| and [IT| present a complete and up to date review, to our knowledge, of the experimental and theoretical data 
on conductivity anisotropy in VO2, including the effect of externally applied stress. Tablejll in cludes the results for 
the conductivity anisotropy, above and below Tmi, when no external stress is appliecP^ISslSniinHsIl^ Table [n| presents 
several results for the strain induced variation ofTj^/ and of the conductivity in VO2 , for situations where hydrostatic 
or uniaxial pressure is applied to the sampleJil^EilillD' 29 40 41 65 67 

As seen in tables |l] and |lTj in most situations where uniaxial pressure is applied the conductivity is measured along 
Cr, the only axis along which it is well defined due to geometry constraints of the samples or to cracking. Few studies 
determine the conductivity along or b^ as a function of applied pressure. In general, it is agreed that (i) an applied 
uniaxial compressive (tensile) stress along the c^-axis leads to an increased (decreased) conductivity at P > Pa//; 
(ii) an applied uniaxial compressive (tensile) stress along the c/^-axis leads to a decreased (increased) Tmi (the axis 
along which Tmi is measured is not always specified); (iii) a small uniaxial stress applied along the [110]/j direction 
has no significant effect on the conductivity nor Tmi along cr but it promotes a phase transition between different 
monoclinic structures (Ml, M2, T/M3, M4), in the insulating phasdSStMl Further experiments are needed in order 
to systematically measure the conductivity along a/? or b/j under (i) an applied uniaxial stress along a/j or b/j; (ii) 
an applied uniaxial stress along c/j; (iii) hydrostatic pressure. Also, clear criteria for distinguishing Tmi from Tst 
would be extremely valuable to help draw a more accurate and complete phase diagram, in line with what has been 
attempted in previous studie^^^^^. The highly oriented strain across a quasi three dimensional structure, achieved in 
epitaxially grown VO2 nanosheets, offers an extraordinary versatility and potential for investigation of these issues. 

Theoretically, several pictures have been suggested to explain the anisotropy in (unstrained) VO2: (1) a two band 
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model description of 3d electrons, within the framework of the Goodenough model, predicts a non conducting ab^j- 
plane for T > Tmi, any residual conductivity in that plane being due to the overlap of 02p and V3d orbitals^^ 
(2) LDA calculations by Allen's group predict the structural distortion to be the main force driving the MIT in 
VO2, following a simple Peierls pictur^^ (3) a three-band Hubbard model, suggested by Tanaka et al, predicts a 
one dimensional conducting phase along c^ for T < Tmi and an isotropically conducting phase for T > Tm^^, (4) 
LDA+cDMFT calculations by Biermann et al. suggest that electron correlations within the e^ levels are weaker than 
those along the aig ones'^; (5) subsequent LDA+cDMFT calculations by Kotliar's group, which include a moderate 
degree of uniaxial strain, suggest that electronic correlations are the main driving force in the MIT, and that the rutile 
phase itself should be able to support both metallic and insulating electronic behavioi!^ (6) Liebsch et al. compared 
the LDA+U, DMFT and GW methods and found that none of them provides a full description of VO2, namely of 
the development of the insulating gap below Tji//™. Nonetheless, most of the calculations based on unstrained VO2 
samples agree that (i) aig and e^ states can be regarded as nearly independent; (ii) the conductivity is expected to 
be of predominantly d-ig behavior for T < Tmi and nearly isotropic for T > Tmi', (hi) electron correlations should 
be included in the model, mainly in the Ci^ oriented aig levels. Our results agree with (i) and (iii), while strain 
induced cracking prevents the accurate description of the effect mentioned in (ii) . Further theoretical investigations 
should take into account not only strairP but also the oxygen degrees of freedom, in order to provide a more accurate 
description of the 02p and V3d orbitals overlap, most relevant for describing the conductivity in the ab^-plane, 
and of the subsequent cracking along cr. Such estimates would be adequate for a quantitative, rather than merely 
qualitative, comparison with our experimental data. 

VI. CONCLUSION 

In summary, we have observed a large anisotropy of the properties of strained lOOnm and 250nm thick VO2 
nanosheets. The increased value of T'^^j compared to Tgt is a clear signature of a Mott- rather than Peierls-driven 
MIT along cr. An orbital tuning picture is proposed to explain the reduced value of the high temperature atj, , 
although a crack induced conductivity decrease cannot be conclusively ruled out. Additional experiments are needed 
in order to clarify the mechanism of the phase transition in strained VO2, mainly in the direction perpendicular to 
Cr. This would allow the development of a more comprehensive phase diagram for this material. In general, epitaxial 
strain engineering is a powerful tool which has the potential to enable careful tuning of the metal-insulator transition 
in numerous other correlated electron materials, thereby providing a viable route towards technologically relevant 
multifunctionality and increased understanding of the microscopic origin of the MIT. 
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FIG. 1. (a) THz transmission along ch and b_R (vertically offset) in the lOOnm (lOO)jj VO2 sample (solid) and in the Ti02 
reference substrate (dashed), below (blue) and above (red) Tmi- The large refractive index anisotropy of the Ti02 substrate 
leads to a slower propagation of the laser pulse along c_r than along b_R, which enables the distinction of the two signals, and 
thus the orientation of the sample with respect to the incident field polarization, (b) THz transmission along c_r and b_R in the 
250nm (lOO)ij VO2 sample, below (blue) and above (red) Tmi, normalized to the high temperature value. The relative low 
temperature transmission along b_R (~ 85%) is dramatically different from that along c_r (~ 15%). 
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FIG. 2. Characterization of the 250nm thick (lOO)j; VO2 film, (a) AFM phase image (0°-5° scale) and corresponding height 
profile: the 250nm sample shows buckling induced ridges along c^ (height indicated by the arrows: Az — 14.307nm). (b) SEM 
image of a section of the sample (seen from the edge, the surface and the cross section being located as labeled) showing a 
~ 30nm wide crack. 
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FIG. 3. Characterization of VO2 thin films on a (lOO)j? Ti02 substrate. Optical images of the (a) 250nm and (b) lOOnm 
films: the samples show buckling induced ridges along c^, spaced by ~ l/im. (c) Temperature dependence (for increasing and 
decreasing temperature) of the a-axis lattice spacing, deduced from XRD data: a ~ 0.1% increase is observed along slr, in the 
lOOnm thick VO2 sample, across the structural transition which occurs at Tst = 340K. The expected linear increase of the 
aij-axis lattice spacing of Ti02 with temperature is also observed. 
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FIG. 4. (a) Room temperature XRD of the lOOnm thick VO2 thin film, along three different directions: (200) ij, (110) ij, (lOl)i;. 
The lattice parameters can be estimated as 4.52A along a.R, 4.46A along b_H and 2.89A along c_r, yielding mismatches of —0.83% 
along a_B, —2.17% along bij and 1.41% along cr. (b) Temperature dependent XRD of the lOOnm thick VO2 thin film. Both 
the Ti02 and the VO2 (200) h peaks are seen to shift as a function of temperature. There is no evidence of the development 
of any additional structural phase. 
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FIG. 5. (a) Temperature dependence of the far infrared conductivity in lOOnm (lOO)ij VO2: — 30(J(,^ above the MIT 
temperature; T^" = 365K while = MOK = tK'=. (b) The VO2 rutile unit cell, following Eyer#i. (c) Effect of tensile 
strain along cr on the electronic structure of metallic VO2: antibonding e^ bands are shifted upwards while the aig band is 
narrowed. 
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FIG. 6. Temperature dependence of the far-infrared conductivity in 250nm (100)_b VO2: (Jcji — 5250(ncm)~^ for T > ~ 
365is:; (76^ < 100{ncmy^ and T^^j can only be estimated as T^Pj ~ SAOK. 
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FIG. 7. Phase diagram of VO2 for stress apphed mostly along [110]. The black solid line indicates Tmi', the white oval defines 
the phase space of our samples; uniaxial stress results are presented in the region delimited by the dashed white lines; doping 
(Vi_^M^02) leading to reduction (M=Nb^+) or oxidation (M=Cr^+) of is specified on the top horizontal axisP^^ 
(refer to the main text for more details). 
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Everhart et al. 


bulk single crystal 
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(1968) 
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Verleur et al. 


bulk single crystal 


reflectivity 
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0.79-2 


(1968) 




(0.25-5eV) 






Continenza et al. 


(calculations) 


model GW 


3-6.7 


N/A 


(1999) 




(O-lQeV) 






Mossanek et al. 


(calculations) 


LDA (0-12eV) 


0.6-> 3 


0.7-3.3 


(2007) 










Lysenko et al. 


30nm film on 


optical (400nm) 


1.05 
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(2007) 


(0121 AhOt. 


diffraction 






Lu et al. 
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5.14 
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(2008) 


(Oil)/? Ti02 
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Tomczak et al. 


(calculations) 
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(2009) 
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lOOnm film on 
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TABLE I. Compilation of previous experimental and theoretical results for the conductiv ity (DC , u nless otherwise specified) 
anisotropy, or in VO2, above and below Tmi, when no external stress is applied.'^^ElHSlIlSHll 
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publication 


sample 
type 


applied 
pressure 


conductivity 
vs. pressure 


conductivity 
along... 


clTMl/dP 

[K/kbar] 


dTMl/dP 
along... 


Minomura et al. 


bulk 


uniaxial 
(Drickamer cell, 
up to 160kbar) 


not specified 


not specified 


-0.46 


not specified 


Neuman et al. 

( 1 QR/l 


single crystal 


hydrostatic 
(up to 6kbar) 


increases 
below Tj\!j 


not specified 


no significant 
change 


not specified 


powder 


hydrostatic 
(up to 6kbar) 


increases 


not specified 


no significant 
change 


not specified 


Berglund et al. 
(1969) 


bulk single crystal 
(0.25x0. 25x0. 2in, 
crack along cr) 


hydrostatic 
(up to 44kbar) 


nicreases ; 
saturates above 
1 Ai I J lo-iUKoar 


Cr 


U.UoZ 


Cr 


Ladd et al. 
(1969) 


bulk single crystal 
iUxiximm ) 


hydrostatic 
(up to 8kbar) 


increases 
(up to oUKDai , 

at Troom) 


CR 


0.06 


CR 


along a J? or b/j 


N/A 


CR 


no significant 
change 


CR 


along CR, 
(up to O.Skbar) 


increases 
b elow Tj\[ J 


CR 


-1.2 


CR 


Pouget et al. 
(1975) 


bulk single crystal 
(0.5mm^x2inni, 
3nim^x4nnn) 


along [nO]R 
(up to 1.2kbar) 


no significant 
change 


CR 


no significant 
cliange 


CR 


Gregg et al. 
(1997) 


thin films on 
Al2O3(012) 


in- plane 
(3-point bend) 


increases 
below Tim 


not specified 


N /A 


N /A 


Muraoka et al. 


10-15nm thick 
single crystal thin 
films on (OOl)fi 
and (110)r TiOa 


along Cr 
(epitaxial! 
-0.3% for {001)r, 
1.2% for (110)/i) 


increases 


aR or bR 
lor [uuijiv. , 
not specified 
for \1W]r 


< 


an or bR 
lor [uuijR, 
not specified 
for [1101 „ 


Arcangeletti et al. 

l^ZUU 1 ) 


single crystal 
(5//m thick slab) 


uniaxial 
(diamond a. c, 
up to 140kbar) 


increases below 

J- MI 

(750-6000cm-^) 


not specified 


N/A 


N/A 


Lu et al. 

l^ZUUo ) 


single crystal 
thin film on 
(Oll)i, TiOj 


in- plane 
(epitaxial! 
-1.2%//[011]fl, 
-0.4%//(0011i,) 


N/A 


an or hn 


< 


ajj or bfl 


Sakai 
(2008) 


polycrystalline thin 
film on metallic Ti 
{aR ill plane) 


in- plane 
(point contact, 
up to 255kbar 


increases 


not specified 


> u 


not specified 


Cao et al. 
(2009) 


single crystal 
nm, along cr) 


along Cr 
(3-pomt bend, 
up to 12kbar) 


increases 
(optical) 


CR 




CR 


our results 
(2010) 


single crystal 
thin film on 
(lOO)iv. Ti02 


in- plane 
(epitaxial! 
5.4%//[001]i?, 
-2.3%//f010lR) 


no significant 
change (far IR)) 


CR 


< 


CR 


decreases above 
Tmi (far IR) 


bfi 


no significant 
change 


bR 



TABLE II. Compilation of previous experimental and theoretical results for the variation of Tmi and of the conductivity (DC, 
unless otherwise specified) in VO2, under applied hydrostatic or uniaxial pressures. We assume AP > for compression and 
AP < for expansion.EI3™ZI^ntelIllM 



